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Abstract 

Using coordinate Bethe ansatze we construct two-magnon states for the sug- 
gested by A. K. Kolezhuk and H.-J. Mikeska family of spin-ladder models with exact 
singlet-rung vacuum. The explicit formula for zero-temperature Raman scattering 
cross section is derived. The corresponding line-shapes are strongly asymmetric and 
their singularities originate only from bound states. This form of a line-shape is in 
a good correspondence to the experimental data. 

1 Introduction 

Raman scattering in spin- ladders was studied in a number of papers (see and 
references therein). The obtained experimental data was analyzed by several theoretical 
approaches [S1,[ZI,[HI- However in none of these papers the exact formula for the Raman 
cross section was used. In the present paper we obtain the exact formula for the special 
class of spin-ladder models with exact singlet-rung vacuum. This family of models was 
first suggested in [TU]. The corresponding Hamiltonian Ti. has the following form: 
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Here Sj,„ (z = 1,2; n = — 00...00) are spin-| operators associated with cites of the ladder 
and / is an identity matrix. The auxihary term H^°n^i in (2) is need only for normalization 
to zero the lowest eigenvalue of the 16 x 16 matrix H of rung-rung interaction. 
It was shown in [TI^ that when the following conditions 
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> 2J||, > -Jc, Jx + J|| > -Jc 



(4) 



are satisfied then the lowest (zero eigenvalue) eigenstate of H is w ^ w, where w is the 
rung-singlet state. In this case the ground state of the Hamiltonian (1) has the simple 
tensor-product form: 

\0) = ll^Wn. (5) 

n 

In order to obtain the full spectrum of H we shall also define the following triplet 
states: 

All other eigenstates of H are separated into the following sectors: singlet ® f^, 
triplet Eijkf^ ® f^, quintet Ujkif^ ® /'^ with eigenvalues: Eq = J± — 2J\\, Si = J± — 
J|| — |Jc, £2 = J± + J\\ — ^Jc, and two triplets w ® f'^ ± ® w with eigenvalues: 
£± = ~ \ Jc ± Jc)- Here tijki = Sik6ji + SuSjk — f^ij^M- 

The Hamiltonian (l)-(3) commutes with the following magnon number operator Q = 
J2n. Qny where Qn = |/+Si,„ -82,71 is the associated with the n-th rung projection operator 
on triplet states. 
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2 The two-magnon states 

Corresponding to (l)-(4) one-magnon states were obtained in jTU]. Suggesting the follow- 
ing Bethe form for two-magnon states \S,f3) (where S is the total spin and (3 the list of 
additional parameters) : 

oo oo 

|0;/3) = ao{m,n;P)...Wm-ifiwm+i---Wn-ifiWn+i..., (7) 



m=— oo n=m+l 
oo oo 



= X] ai{m,n; P)eijk ■ ■ .Wm-ifmWm+i ■ ■ -Wn-ifnWn+i . . . , (8) 

m=— oo n=m+l 
oo oo 

|2;/3)ji = a2{m,n; P)tijkl . . .Wrn-lf!^Wrn+l ■ ■ -Wn-lfnWn+l . ■ . (9) 

m=— oo n=m+l 

we obtain in standard way ^T] the following Schrodinger equation: 

^[as{m - 1, n; (3) + as{m + 1, n; (3) + as{m, n - 1; P) + as{m, n + 1; 

+ {2J^ - 3J,)as{m,n; (3) = Eas{m,n; (3), (10) 
and Bethe condition for amplitudes: 

2Asasin, n + 1; p) = as{n, n; (3) + as{n + 1, n + 1; /?). (11) 
Here A— ^^~^+~^ 



£+ — £_ 

For each S the Eq. (11) has two solutions. The scattering solution: 

ar"(m, n- h, k2) = C5,i2e'(''i'"+'^^") - C7s,2ie'('=^™+'=^"), (12) 
with Cs,ab = COS '""^ '' — As'C* 2 , and the bound solution: 

where the real parameters v > and —it < u < n satisfy the following condition: 

cos^t = Ase"^ (14) 
From (14) and nonnegativity of v follows that 

|cosw| < |As| <e''. (15) 
The corresponding to (12) and (13) eigenvalues are: 

^|™**(A;i,fc2) = 2Jx-3Jc + Jc(cos/ci + cosA;2), (16) 

^feoW(^) = 2Jx+ (As-3)Je+ ^cos^w. (17) 

A5 

As we see from (12) and (13) the translation invariant states correspond to a|'^'^**(m, n; k, —k), 
a^°™'^(m, n; 0) and a^°""'^(m, n; tt). 



3 Calculation of Raman cross section 

Following Sugai we shall consider only the case when the incident and scattered light 
have parallel polarization directions both lying in the plane of the ladder and forming an 
angle 9 with respect to vertical bonds (rungs). The zero-temperature two-magnon Raman 
scattering cross section as a function of frequency and 9 can be expressed using Fermi's 
golden rule 

J(..,^) = Jm^^^5^|(«|7^«(^^)|0)r<5(^-E„), (18) 

a 

where 2N + 1 is the number of rungs. Within the Fleury-Loudon-Elliot approach the 
effective Raman Hamiltonian 7i^{9) have the following form pP,[Sj (we also take into 
account interactions across diagonals): 

H''{9) = Ai,gCos^9n''' + Ad.ag{cos\e+-f)n''' + cos\e--f)n''') + ArungSm^9n'^^'. (19) 

Here Ai^g, A^iag and Arung are constants and 7 is the angle between rung and diagonal 
directions. Operators H^'"'^ and are the following: 

'j_^rung ^ ^ g^^^ _ g^^^^ 'j_^leg ^ ^ g,^^ _ Sj,ri+1, TY'^^'^' = ^ Si(2),n " 82(1),™+!. (20) 

n i,n n 

Expressing 7i'^^, Ti'^^, Ti'^'^ and 7-^^""9 from the auxiliary operators: 

'H'^^ = ^(Si^n ± 82,™) ■ (Si^„+i ± S2,n+l). (21) 

n 

and taking into account the Eq. (6) we represent I{uj, 0) in the factorized form: 

/(W, 9) = ^Mleg + Adiag Sm^ 7) sin^ d + Adiag COS^ 7 COS^ 9f Io{uj) , (22) 

where 

/o(a;)= lim |(a|H--|0)p<5(^ - i?„). (23) 

a 

Formula (22) expresses the polarization angle dependence of Raman cross section 
however it may be applied in a straightforward way only for 9 = 0. 

From the Eq. (6), translational and SU(2) invariance of Ti. follows that only trans- 
lation invariant singlet two-magnon states contribute to the formula (23). Separating the 
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contributions from scattering and bound states we obtain: 

ir'\u;) = lim y \^n=-N<nn + i-k,~k)\' s{uj - Er'\k, -k)), (24) 

j6oW(^) = lim y l^n=-iv"(^'^ + ^'^)l' (25) 
From (12) and (13) follows: 

N n-1 

Yl Yl K'''\rn,n;k,-k)\^ = 4N\l - 2 Ao cos k + Al) + 0{N), (26) 

m=-N 
N 

J2 a^'"**(n,n + l;A;,-fc)| = 4Nsmk + 0{l), (27) 

" ' m=-N 
N 

a^'^'^'^{n,n + l-u)\ = {2N + l)e-\ u = 0,7r. (29) 



n=-N+l m=-N 
N 



n=-N+l m=-N 
N 



n=-N 



Using the substitution Y^j^ lo^ ^k we obtain from (26)-(29) the final expressions 

for the cross sections: 



'° - Ul + AI- 2.A0) ' 

/o^-V) = ^(l-^)e(A^-l)5(2x-Ao-i-). (31) 

Here © is the step function and x — '^"^27"'"^'^° rescaling parameter. 

Prom (15) and (31) follows that the contribution from bound states exist 
only for | Ao| > 1. The behavior of /^"^"** as a function of x also essentially depends on the 
parameter Aq = ^ — 2-^^. When Aq = ±1 the formula (30) reduces and the hne-shape 
has a singularity at x = Aq. For Aq = 1 it lies in the top of the two-magnon continuum 
however for Aq = — 1 in the bottom. For Aq 7^ ±1 the cross section J*'^"** is a regular 
function of x and has the maximum in the point Xmax = -t^zt- 

In order to study the line-shape in more detail we shall find its inflection points. Cal- 
culating the second derivative of J^'^"** with respect to x we obtain the following condition: 

p{x, Ao) = 4Ao(l + A^)^^ - 12A^x' - A^ + GA^ - 1 = 0. (32) 
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Since p(±l, Aq) = —(1 ^AqY, the polynomial p{x, Aq) for Aq 7^ ±1 has only or 2 zeros 
in the interval (—1, 1). From standard calculation follows that p{x, Aq) has the maximum 
Pmax = — Aq + 6Aq — 1 in the point x = 0. It is evident now that for Pmax > the 
line-shape of J*™** has two inflection points. From the straightforward calculation follows 
that Pmax > only for 

A_ < |Ao| < A+, (33) 

where A± = \/3 ± 2^2 (A_ ^ 0.4142, A+ ^ 2.4142). It may be easily proved in a 
straightforward way that (A+ — A_)^ = 4, so A+ — A_ = 2. 

In the case (33) the line-shape near the Xmax is similar to van-Hove singularity. For 
A_ < Aq < A+ this "singularity" lies near the top of the two-magnon continuum however 
for — A+ < Aq < — A_ near the bottom. In both the cases the line-shape of Raman 
scattering is strongly asymmetric. The case Pmax < with no inflection points may be 
interpreted as a broad maximum. Some line shapes corresponding to different values of 
Aq are presented in the Fig. 1. 

As it follows from (16) and (17) for Aq —>■ ±1 + 0^ the top (bottom) of the two-magnon 
continuum and the bound two-magnon state have the same energy: 2 J± — 3 Jc ± 2 Jc. It 
was proposed in [Hj that in this case the resonance between bound and scattering states 
leads to a redistribution of Raman intensity and merging of singularity. However as we 
see from (30) and (31) in our model this conjecture fails. Moreover the singularity in /^™** 
appears only in the resonance Aq = ±1 cases. 

4 Comparison with experiment and discussion 

Raman scattering in MgVgOs and CaV205 were reported in P^. It was pointed that 
for both materials the corresponding line-shapes are strongly asymmetric and have one 
maximum instead of two. This fact was considered as strange and there were suggested 
some conjectures to interpret it. For example it was supposed in that in MgVaOs 
there is no spin-gap and the magnetic ordering is 2D, or the spin-gap is so small (about 
10 cm~^) that can not be observed by the used experimental resolution. The asymmetry of 
the line-shape for CaV205 was interpreted in jH] as originating from next-nearest neighbor 
interactions. In [Zj it was conjectured that the second peak in the line-shape of CaV205 
is not observed because it is dominated by phonon peak. In it was conjectured that 
the asymmetry originates from resonance with two-triplet bound state. 
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In our paper we have demonstrated that the hne-shape asymmetry in spin-ladder Ra- 
man scattering is not something strange and outstanding but may appear in a sufficiently 
big class of models. Of course we do not pretend that for some values of exchange param- 
eters our toy model necessary describes the real materials such as CaV205 or MgVaOs. 
Nevertheless perhaps the true ground state is in some sense "close" to our idealized one 
(5) and we may believe that our model correctly represents some general qualitative fea- 
tures of real materials. In this context we emphasize that the exactly calculated Raman 
scattering line-shape may be strongly asymmetric without any additional assumptions 
such as next-nearest neighbor interactions, resonance with bound state or dominating by 
phonon peak. 

The author is very grateful for discussions with K. P. Schmidt, P. P. Kulish, and 
L. V. Prokhorov. 
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